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Summary. — We introduce a method to measure with high resolution the solar 
diameter from the ground, through the eclipse observations by reconsidering the 
definition of the solar edge. The outer part of the Limb Darkening Function (LDF) 
is recovered using the luminosity evolution of a Baily's Bead and the profile of the 
lunar limb available from the Kaguya satellite. The method proposed is applied for 
the videos of the eclipse in January, 15, 2010 recorded by Richard Nugent in Uganda 
and Andreas Tegtmeier in India. The result shows light from solar limb detected 
at least 0.65 arcsec beyond the LDF inflection point, and this fact may suggest to 
reconsider the evaluations of the historical eclipses made with naked eye. 

PACS95.10.Gi - Eclipses. 96. 60.-j - Solar physics. 96.60.Bn - Solar diameter. 



1. The method of eclipses 

With the echpse observations we are able to bypass some problems that affect the 
measurement of solar diameter. The atmospheric and instrumental effects that distort 
the shape of the limb [1] are overcome by the fact that the scattering of the Sun's light 
is greatly reduced by the occultation of the Moon, therefore there are much less photons 
from the photosphere to be poured, by the PSF effect, in the outer region. 

I'l. Comparison between observation and ephemeris. - The method exploits the ob- 
servation of the beads of light that appear or disappear from the bottom of a lunar valley 
when the solar limb is almost tangent to the lunar limb. 

The timing of their appearing or disappearing and their positions are measured. These 
instants when the photosphere disappears or emerges behind the valleys of the lunar 
limb, are determined solely by the positions of the Sun and the Moon with respect to 
the observer; their angular size depends on the lunar marginal profile while the effects of 
the atmospheric seeing on these phenomena are negligible. 
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Bead's magnitude evolution 




0.05 



0.1 



0.15 



0.2 



Fig. 1. - Beads magnitude evolution: the height of the solar limb above the valley is on the 
abscissa. The various square dots represent different type of telescopes. [25, p] means 25 cm 
opening with projection of the image. [25, ND4] the 25 cm telescope with a filter of Neutral 
Density of transmittance 1/10000; ND5 stands for 1/100000, and so on. 



The International Occultation Timing Association (IOTA) is currently engaged to 
observe the eclipses with the aim of measuring the solar diameter. This is facilitated by 
the development of the software Occult 4 by David Herald(^). 

The technique consists to look at the time of appearance of the beads and to compare it 
with the calculated positions by the ephemeris using the software Occult 4. The simulated 
Sun by Occult 4 has the standard radius: 959.63 arcsec at lAU [2]. The difference 
between the simulations and the observations is a measure of the radius correction with 
respect to the standard radius (AR). 

In the lunar polar regions, due to the geometry of the eclipse, the beads last for more 
time [3] . Observing in the lunar polar regions we are also able to avoid the measurement 
in the solar active regions that could affect the measures [4]. In fact the solar active 
regions do not appear in solar latitude higher than ^ 40°. The maximum offset between 
the lunar and the solar pole is ^9°. 

However this approach conceives the bead as an on-off signal. In this way one as- 
sumes the Limb Darkening Function as a step function, but it is actually not. Different 
compositions of optical instruments (telescope + filter + detector) could have different 
sensitivities and different signal/noise ratios, recording the first signal of the bead in cor- 
respondence of different points along the luminosity profile. This leads to different values 
of the perceived AR i.e. the correction to the standard angular solar radius at 1 AU (see 
Fig. 1). An improvement of this approach has to take into account the whole shape of 
the Limb Darkening Function and thus the actual position of the inflection point. 

1'2. Numerical calculation. - The shape of the light curve of the bead is determined 
by the shape of the LDF (not affected by seeing) and the shape of the lunar valley that 
generates the bead. Calling w(x) the width of the lunar valley (i.e. the length of the solar 
edge visible from the valley in function of the height x from the bottom of the valley) , 
and B(x) the surface brightness profile (i.e. the LDF), one could see the light curve L(y) 



(^) www.lunar-occultations.com/iota/occult4.htm 
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as a convolution of B(x) and w(x), being | y \ the distance between the botton of the 
hmar valley and the standard solar edge, setting to the position of the standard edge. 
L{y) = J B{x) w{y — x) dx 
The discrete convolution is: 

L{m) = ^ B[n) w{m — n)h 

where n, m are the index of the discrete layers corresponding to x, y coordinate and h 
the layers thickness. 

Thus the profile of the LDF is discretized in order to obtain the solar layers B(n) of 
equal thickness and concentric to the center of the Sun. In the short space of a lunar 
valley these layers arc roughly parallel and straight. 

The lunar valley is also divided in layers of equal angular thickness. During a bead 
event every lunar layer is filled by one solar layer every given interval of time. Step by 
step during an emerging bead event, a deeper layer of the solar atmosphere enters into 
the profile drawn by the lunar valleys (see Fig. 2), and each layer casts light through the 
same geometrical area of the previous one. 

Being: Ai, A2.. An the area of lunar layers, from the bottom of the valley going 
outward; Bi, B2-. Bn the surface brightness of the solar layers (our goal) from the outer 
(dimmer) going inward; ii, L2.. the value of the observed light curve from the first 
signal to the saturation of the detector or to the replenishment of the lunar valley. One 
has: 

Lm = BiA„i + B2Am-l + .. + B„iAi = 

To derive the LDF profile: 
Bi = iiMi 

B2 = [L2 - {B, ■ A2)] Ml 

B3 = [L3 - {Bi ■ As + B2 ■ A2)] Ml 

Bi - [L4 - (Bi ■A4 + B2-A3 + B3- A2)] Ml 

and so on. 

The situation described above relates to an emerging bead. The same process can run 
for a disappearing bead, simply plotting the light curve back in time. 
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2. An application of the method of eclipses 

We studied the videos of the annular eclipse in January 15, 2010 realized by Richard 
Nugent in Uganda and Andreas Tegtnieier in India. 

The equipment of Nugent was: CCD camera Watec 902H Uhimate(^); Matsukov 
telescope (90/1300 mm); panchromatic filter Thousand Oaks, ND5. 

The equipment of Tegtmeier was: CCD camera Watec 120N(3); Matsukov telescope 
(100/1000 mm); Fiher: I = lOTA/ES green glass based neutral ND4. 

2'1. Analysis. - Two beads located at Axis Angle(^) (AA) — 171°andl77° are analyzed 
for both the videos. 

The lunar valley analysis is performed with the new lunar profile obtained by the 
laser altimeter (LALT) onboard the Japanese lunar explorer Kaguya(^). The radial 
topographic error is estimated to be 4.1 m corresponding to 2.1-2.3 mas at the distance 
of the Earth. The choice of the thickness of the layers has to be optimal: large enough 
to reduce B„ uncertainties, but small enough to have a good resolution of the LDF. We 
chose h = 30 mas for the lunar valley at A A = 177° and h ^ 73 mas for the lunar valley 
at AA = 171°. 

The motion of the Sun in the lunar valley is simulated with the Occult 4 software. 
A program in Fortran90 is performed to calculate the LDF points (i?n), taking into 
account the light curves of the beads, the lunar profiles and the respective uncertainties. 

2'2. Results. - The resulting points show that the inflection point is clearly between 
the two profile obtained for each of the two beads. The saturation of the CCD pixels 
avoided to measure the luminosity function more inward for Nugent's video, while the 
low sensitivity avoided to measure the luminosity function more outward for Tegtmeier's 
video. Therefore from these points it is impossible to infer an exact location of the 
inflection point. But it is possible to deduce an upper and lower limit corresponding to 
the points that better constrain it: -0.190 arcsec < Ai? < -1-0.050 arcsec. 

The analysis of these beads with the "classical approach" would lead to a difference 
greater than one arcsecond for AR between the 2 videos. The improvement introduced 
by this approach is then shown. 

3. Historical eclipses 

Even if observed with the naked eyes, historical eclipses are of great interest for 
studying the variations of the solar diameter over periods longer than a solar cycle. 

As an example in April 9, 1567 in Rome, Chrispopher Clavius observed an annular 
eclipse [6]. The solar limb should have been at least Ai? > -1-2.5 arcsec in order to be 
higher than the mountains of the Moon. Eddy et al.[7] deduced from this interpretation 
a secular shrinking of the Sun from the 16th century to the present. In general the 
evaluations of historical eclipses is made without any reference to the inflection point, 
and considering the LDF as a Heaviside profile. This can enlarge the measured AR. It 
is the main concern with naked eye observations. 



(^) http: / /www. aegis- elec.com/products/watec- 902H_spec_eng.pdf 

(^) http:/ /www. aegis-elec.com/products/watec- 100N_spec_eng.pdf 

(*) the angle around the limb of the Moon, measured Eastward from the Moon's North pole 

(^) http: / /wms. selene.jaxa.jp/selene_viewer/index_e. html 
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Fig. 3. - The luminosity profiles obtained for the bead at AA = 177° are plotted and put 
together. The inner and brighter part is obtained from Tegtmeier's video; the outer and weaker 
part is obtained from Nugent's video. The luminosity profile is normalized to the center of the 
solar disk according to Rogerson [5] for the inner part, and in arbitrary way for the outer part. 
The zero of the abscissa is the position of the standard solar limb with a radius of 959.63 arcsec 
at 1 AU. The error bars on y axis are the 90% confidence level. The error bars on x axis are 
the thickness (h) of the lunar layers. The solid line is an interpolation between the profiles and 
gives a possible scenario on the position of the infiection point. 



4. Conclusions 

Wc obtained a detailed profile of the outer part of the LDF, demonstrating tlie func- 
tionality of the method. Although it was impossible to observe the inflection point, a 
consideration is coming out: if the LDF is considered as a step function, like up to now 
in the eclipse analyses, the solar diameter can be overstimated. And this effect can apply 
to historical eclipses, of which we reconsidered the Clavius' one (1567). Moreover the 
LDF profile, around the inflection point position, is an important source of information 
for the solar atmosphere [4]. Improvements to obtain a detailed profile comparable to 
the space astrometry arc possible: 

• An increased dynamic range of the CCD detectors from 8 bits to 12 bits is recom- 
mended in order to avoid CCD saturation before reaching the inflection point. 

• The profile is dependent on the wavelength [1,4]. One must then observe at specific 
photometric bands to be able to compare the results with other measures. 
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